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(57) Abstract: Techniques to calibrate the downlink and uplink channels to account for differences in the frequency responses of 
the transmit and receive chains at an access point and a user terminal. In one embodiment, pilots are transmitted on the downlink 
^5 and uplink channels and used to derive estimates of the downlink and uplink channel responses, respectively. Two sets of correction 
2 factors are then determined based on the estimates of the downlink and uplink channel responses. A calibrated downlink channel 
is formed by using a first set of correction factors for the downlink channel, and a calibrated uplink channel is formed by using a 
second set of correction factors for the uplink channel. The first and second sets of correction factors may be determined using a 
matrix-ratio computation or a minimum mean square error (MMSE) computation. The calibration may be performed in real-time 
based on over-the-aii 
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CHANNEL CALIBRATION FOR A TIME DIVISION DUPLEXED 
COMMUNICATION SYSTEM 

Claim of Priority under 35 U.S.C. §119 
[0001] This application claims the benefit of U.S. Provisional Application Serial No. 

60/421,462, entitled, "Channel Calibration for a Time Division Duplexed 
Communication System," and U.S. Provisional Application Serial No. 60/421,309, 
entitled, "MIMO WLAN System," both of which are filed on October 25, 2002, 
assigned to the assignee of the present application, and incorporated herein by reference. 

BACKGROUND 

Field 

[0002] The present invention relates generally to communication, and more specifically 

to techniques for calibrating downlink and uplink channel responses in a time division 
duplexed (TDD) communication system. 

Background 

[0003] In a wireless communication system, data transmission between an access point 

and a user terminal occurs over a wireless channel. Depending on the system design, 
the same or different frequency bands may be used for the downlink and uplink. The 
downlink (or forward link) refers to transmission from the access point to the user 
terminal, and the uplink (or reverse link) refers to transmission from the user terminal to 
the access point. If two frequency bands are available, then the downlink and uplink 
may be transmitted on separate frequency bands using frequency division duplexing 
(FDD). If only one frequency band is available, then the downlink and uplink may 
share the same frequency band using time division duplexing (TDD). 

[0004] To achieve high performance, it is often necessary to know the frequency 

response of the wireless channel. For example, the response of the downlink channel 
may be needed by the access point to perform spatial processing (described below) for 
downlink data transmission to the user terminal. The downlink channel response may 
be estimated by the user terminal based on a pilot transmitted by the access point. The 
user terminal may then send the channel estimate back to the access point for its use. 
For this channel estimation scheme, a pilot needs to be transmitted on the downlink and 
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additional delays and resources are incurred to send the channel estimate back to the 
access point. 

[0005] For a TDD system with a shared frequency band, the downlink and uplink 

channel responses may be assumed to be reciprocal of one another. That is, if H 
represents a channel response matrix from antenna array A to antenna array B, then a 
reciprocal channel implies that the coupling from array B to array A is given by H r , 
where M r denotes the transpose of matrix M . Thus, for a TDD system, the channel 
response for one link may be estimated based on a pilot sent on the other link. For 
example, the uplink channel response may be estimated based on an uplink pilot, and 
the transpose of the uplink channel response estimate may be used as an estimate of the 
downlink channel response. 

[0006] However, the frequency responses of the transmit and receive chains at the 

access point are typically different from the frequency responses of the transmit and 
receive chains at the user terminal. In particular, the frequency responses of the 
transmit/receive chains used for uplink transmission may be different from the 
frequency responses of the transmit/receive chains used for downlink transmission. The 
"effective" downlink channel response (i.e., including the transmit/receive chains) 
would then be different from the reciprocal of the effective uplink channel response due 
to differences in the transmit/receive chains (i.e., the effective channel responses are not 
reciprocal). If the reciprocal of the channel response estimate obtained for one link is 
used for spatial processing on the other link, then any difference in the frequency 
responses of the transmit/receive chains would represent error that, if not determined 
and accounted for, may degrade performance. 

[0007] There is, therefore, a need in the art for techniques to calibrate the downlink and 

uplink channels in a TDD communication system. 

SUMMARY 

[0008] Techniques are provided herein to calibrate the downlink and uplink channels to 

account for differences in the frequency responses of the transmit and receive chains at 
the access point and user terminal. After calibration, an estimate of the channel 
response obtained for one link may be used to obtain an estimate of the channel 
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response for the other link. This can then simplify the channel estimation and spatial 
processing. 

[0009] In one embodiment, a method is provided for calibrating the downlink and 

uplink channels in a wireless TDD multiple-input multiple-output (MTMO) 
communication system. In accordance with the method, a pilot is transmitted on the 
uplink channel and used to derive an estimate of the uplink channel response. A pilot is 
also transmitted on the downlink channel and used to derive an estimate of the downlink 
channel response. Two sets of correction factors are then determined based on the 
estimates of the downlink and uplink channel responses. A calibrated downlink channel 
is formed by using a first set of correction factors for the downlink channel, and a 
calibrated uplink channel is formed by using a second set of correction factors for the 
uplink channel. The appropriate correction factors will be used at the respective 
transmitter for the downlink and uplink channels. The responses of the calibrated 
downlink and uplink channels are approximately reciprocal due to the two sets of 
correction factors. The first and second sets of correction factors may be determined 
using a matrix-ratio computation or a minimum mean square error (MMSE) 
computation, as described below. 

[0010] The calibration may be performed in real-time based on over-the-air 

transmission. Each user terminal in the system may derive the second set of correction 
factors for its own use. The first set of correction factors for the access point may be 
derived by multiple user terminals. For an orthogonal frequency division multiplexing 
(OFDM) system, the calibration may be performed for a first set of subbands to obtain 
two sets of correction factors for each subband in the set. Correction factors for other 
"uncalibrated" subbands may be interpolated based on the correction factors obtained 
for the "calibrated" subbands. 

[0011] Various aspects and embodiments of the invention are described in further detail 

below. 



BRIEF DESCRIPTION OF THE DRAWINGS 
[0012] The features, nature, and advantages of the present invention will become more 

apparent from the detailed description set forth below when taken in conjunction with 
the drawings in which like reference characters identify correspondingly throughout and 
wherein: 
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[0013] FIG. 1 shows the transmit and receive chains at an access point and a user 

terminal in a MIMO system; 
[0014] FIG. 2 illustrates the application of correction factors to account for differences 

in the transmit/receive chains at the access point and user terminal; 
[0015] FIG. 3 shows a process for calibrating the downlink and uplink channel 

responses in a TDD MIMO-OFDM system; 
[0016] FIG. 4 shows a process for deriving estimates of the correction vectors from the 

downlink and uplink channel response estimates; 
[0017] FIG. 5 is a block diagram of the access point and user terminal; and 

[0018] FIG. 6 is a block diagram of a TX spatial processor. 

DETAILED DESCRIPTION 
[0019] The calibration techniques described herein may be used for various wireless 

communication systems. Moreover, these techniques may be used for single-input 
single-output (SISO) systems, multiple-input single-output (MISO) systems, single- 
input multiple-output (SIMO) systems, and multiple-input multiple-output (MEMO) 
systems. 

[0020] A MEVIO system employs multiple (Nt) transmit antennas and multiple (N R ) 

receive antennas for data transmission. A MEVIO channel formed by the N T transmit 
and Nr receive antennas may be decomposed into N s independent channels, with 
N s < min{A^ r , N R } . Each of the N s independent channels is also referred to as a spatial 
subchannel or an eigenmode of the MEVIO channel and corresponds to a dimension. 
The MEVIO system can provide improved performance (e.g., increased transmission 
capacity) if the additional dimensionalities created by the multiple transmit and receive 
antennas are utilized. This typically requires an accurate estimate of the channel 
response between the transmitter and receiver. 

[0021] FIG. 1 shows a block diagram of the transmit and receive chains at an access 

point 102 and a user terminal 104 in a MEVIO system. For this system, the downlink 
and uplink share the same frequency band in a time division duplexed manner. 

[0022] For the downlink, at access point 102, symbols (denoted by a "transmit" vector 

*dn ) ^ processed by a transmit chain (TMTR) 114 and transmitted from N ap antennas 
1 16 over a wireless channel. At user terminal 104, the downlink signals are received by 
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N ul antennas 152 and processed by a receive chain (RCVR) 154 to provide received 
symbols (denoted by a "receive" vector r^). The processing by transmit chain 114 
typically includes digital-to-analog conversion, amplification, filtering, frequency 
upconversion, and so on. The processing by receive chain 154 typically includes 
frequency downconversion, amplification, filtering, analog-to-digital conversion, and so 
on. 

[0023] For the uplink, at user terminal 104, symbols (denoted by transmit vector x up ) 

are processed by a transmit chain 164 and transmitted from N ul antennas 152 over the 
wireless channel. At access point 102, the uplink signals are received by N ap antennas 
116 and processed by a receive chain 124 to provide received symbols (denoted by 
receive vector r up ). 

[0024] For the downlink, the receive vector at the user terminal may be expressed as: 

r dn =R ut HT ap x dn , Eq(l) 
where x^ is the transmit vector with N ap entries for the symbols transmitted from the 
N ap antennas at the access point; 
r dn is the receive vector with N ut entries for the symbols received on the N ut 

antennas at the user terminal; 
T ap is an N ap xN ap diagonal matrix with entries for the complex gains associated 

with the transmit chain for the N ap antennas at the access point; 
R ut is an N ut x N ut diagonal matrix with entries for the complex gains associated 

with the receive chain for the N ut antennas at the user terminal; and 
H is an N m x N ap channel response matrix for the downlink. 
The responses of the transmit/receive chains and the wireless channel are typically a 
function of frequency. For simplicity, a flat-fading channel (i.e., with a flat frequency 
response) is assumed. 

[0025] For the uplink, the receive vector at the access point may be expressed as: 

r up =R ap H r T ut x up , Eq(2) 
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where x up is the transmit vector for the symbols transmitted from the N ut antennas at 
the user terminal; 

r up is the receive vector for the symbols received on the N ap antennas at the 
access point; 

T ut is an N ut X N u[ diagonal matrix with entries for the complex gains 
associated with the transmit chain for the N ul antennas at the user 
terminal; 

R ap is an N ap x N ap diagonal matrix with entries for the complex gains 
associated with the receive chain for the N ap antennas at the access point; 
and 

H T is an N ap x N ul channel response matrix for the uplink. 

[0026] For a TDD system, since the downlink and uplink share the same frequency 

band, a high degree of correlation normally exists between the downlink and uplink 
channel responses. Thus, the downlink and uplink channel response matrices may be 
assumed to be reciprocal (i.e., transposes) of each other and denoted H and H r , 
respectively, as shown in equations (1) and (2). However, the responses of the 
transmit/receive chains at the access point are typically not equal to the responses of the 
transmit/receive chains at the user terminal. The differences then result in the following 
inequality R ap H r T ut * (R ut HT ap ) r . 

[0027] From equations (1) and (2), the "effective" downlink and uplink channel 

responses, and H up , which include the responses of the applicable transmit and 
receive chains, may be expressed as: 

S dn =R ut HT ap and H up = R ap H r T m . Eq(3) 
Combining the two equations in equation set (3), the following relationship may be 
obtained: 

SiHdolJ = (RapHupTu,) 7 " =i; 1 t H;R; 1 p . Eq (4) 

Rearranging equation (4), the following is obtained: 
Hl p =l ut R- u 1 t H dn TjR ap =K; 1 H dn K ap 

or 
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H up =(K: t 1 H dn K ap ) r , Eq(5) 
where K ut = T;' t R ut and K ap = T ap R ap . Equation (5) may also be expressed as: 

H up K ut =(H dn K ap ) r . Eq(6) 
[0028] The left-hand side of equation (6) represents the calibrated channel response on 

the uplink, and the right-hand side represents the transpose of the calibrated channel 
response on the downlink. The application of the diagonal matrices, K ut and K ap , to 
the effective downlink and uplink channel responses, as shown in equation (6), allows 
the calibrated channel responses for the downlink and uplink to be expressed as 
transposes of each other. The (N ap xN ap ) diagonal matrix K ap for the access point is 
the ratio of the receive chain response R ap to the transmit chain response T ap (i.e., 
Kap =7[r L )> where the ratio is taken element-by-element. Similarly, the (N ut xiV ur ) 

■lap 

diagonal matrix K ut for the user terminal is the ratio of the receive chain response R ut 
to the transmit chain response T ut . 

[0029] The matrices K ap and K ut include values that can account for differences in the 

transmit/receive chains at the access point and user terminal. This would then allow the 
channel response for one link to be expressed by the channel response for the other link, 
as shown in equation (6). 

[0030] Calibration may be performed to determine the matrices K ap and K ut . 

Typically, the true channel response H and the transmit/receive chain responses are not 
known nor can they be exactly or easily ascertained. Instead, the effective downlink 
and uplink channel responses, H dn and H up , may be estimated based on pilots sent on 
the downlink and uplink, respectively, as described below. Estimates of the matrices 
K ap and K ut , which are referred to as correction matrices K ap and K ut , may then be 
derived based on the downlink and uplink channel response estimates, H dn and H up , as 
described below. The matrices K ap and K ut include correction factors that can account 
for differences in the transmit/receive chains at the access point and user terminal. 
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[0031] FIG. 2 illustrates the application of the correction matrices K ap and K ut to 

account for differences in the transmit/receive chains at the access point and user 
terminal. On the downlink, the transmit vector is first multiplied with the matrix 
K ap by a unit 112. The subsequent processing by transmit chain 114 and receive chain 
154 for the downlink is the same as shown in FIG. 1. Similarly, on the uplink, the 
transmit vector x up is first multiplied with the matrix K ut by a unit 162. Again, the 
subsequent processing by transmit chain 164 and receive chain 124 for the uplink is the 
same as shown in FIG. 1. 

[0032] The "calibrated" downlink and uplink channel responses observed by the user 

terminal and access point, respectively, may then be expressed as: 

H cdn =H dn K ap and H cup =H up K ut , Eq (7) 

where and H cup are estimates of the "true" calibrated channel response 

expressions in equation (6). Combining the two equations in equation set (7) using the 
expression in equation (6), it can be shown that H cup « H ^ . The accuracy of the 
relationship H cup = is dependent on the accuracy of the matrices K ap and K ut , 
which in turn is typically dependent on the quality of the downlink and uplink channel 
response estimates, and H up . 

[0033] As shown above, calibration may be performed in a TDD system to determine 

the differences in the responses of the transmit/receive chains at the access point and 
user terminal, and to account for the differences. Once the transmit/receive chains have 
been calibrated, a calibrated channel response estimate obtained for one link (e.g., H cdn ) 
may be used to determine an estimate of the calibrated channel response for the other 
link (e.g., H cup ). 

[0034] The calibration techniques described herein may also be used for wireless 

communication systems that employ OFDM. OFDM effectively partitions the overall 
system bandwidth into a number of (N F ) orthogonal subbands, which are also referred to 
as frequency bins or subchannels. With OFDM, each subband is associated with a 
respective subcarrier upon which data may be modulated. For a MIMO system that 



WO 2004/039022 



PCT7US2003/034515 



utilizes OFDM (i.e., a MMO-OFDM system), each subband of each eigenmode may be 
viewed as an independent transmission channel. 

[0035] The calibration may be performed in various manners. For clarity, a specific 

calibration scheme is described below for a TDD MTMO-OFDM system. For this 
system, each subband of the wireless link may be assumed to be reciprocal. 

[0036] FIG. 3 is a flow diagram of an embodiment of a process 300 for calibrating the 

downlink and uplink channel responses in the TDD MTMO-OFDM system. Initially, 
the user terminal acquires the timing and frequency of the access point using acquisition 
procedures defined for the system (step 310). The user terminal may then send a 
message to initiate calibration with the access point, or the calibration may be initiated 
by the access point. The calibration may be performed in parallel with 
registration/authentication of the user terminal by the access point (e.g., during call 
setup) and may also be performed whenever warranted. 

[0037] The calibration may be performed for all subbands that may be used for data 

transmission (which are referred to as the "data" subbands). Subbands not used for data 
transmission (e.g., guard subbands) typically do not need to be calibrated. However, 
since the frequency responses of the transmit/receive chains at the access point and user 
terminal are typically flat over most of the band of interest, and since adjacent subbands 
are likely to be correlated, the calibration may be performed for only a subset of the data 
subbands. If fewer than all data subbands are calibrated, then the subbands to be 
calibrated (which are referred to as the "designated" subbands) may be signaled to the 
access point (e.g., in the message sent to initiate the calibration). 

[0038] For the calibration, the user terminal transmits a MTMO pilot on the designated 

subbands to the access point (step 312). The generation of the MEMO pilot is described 
in detail below. The duration of the uplink MIMO pilot transmission may be dependent 
on the number of designated subbands. For example, 8 OFDM symbols may be 
sufficient if calibration is performed for four subbands, and more (e.g., 20) OFDM 
symbols may be needed for more subbands. The total transmit power is typically fixed, 
so if the MEMO pilot is transmitted on a small number of subbands, then higher amounts 
of transmit power may be used for each of these subbands and the SNR for each 
subband is high. Conversely, if the MEMO pilot is transmitted on a large number of 
subbands then smaller amounts of transmit power may be used for each subband and the 
SNR for each subband is worse. If the SNR of each subband is not sufficiently high, 
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then more OFDM symbols may be sent for the MIMO pilot and integrated at the 
receiver to obtain a higher overall SNR for the subband. 

[0039] The access point receives the uplink MIMO pilot and derives an estimate of the 

uplink channel response, H up (fc), for each of the designated subbands, where k 
represents the subband index. Channel estimation based on the MIMO pilot is 
described below. The uplink channel response estimates are quantized and sent to the 
user terminal (step 314). The entries in each matrix H up (£) are complex channel gains 
between the N ut transmit and N ap receive antennas for the uplink for the fc-th subband. 
The channel gains for all matrices may be scaled by a particular scaling factor, which is 
common across all designated subbands, to obtain the desired dynamic range. For 
example, the channel gains in each matrix H up (£) may be inversely scaled by the 
largest channel gain for all matrices H up (fc) for the designated subbands so that the 
largest channel gain is one in magnitude. Since the goal of the calibration is to 
normalize the gain/phase difference between the downlink and uplink, the absolute 
channel gains are not important. If 12-bit complex values (i.e., with 12-bit inphase (I) 
and 12-bit quadrature (Q) components) are used for the channel gains, then the 
downlink channel response estimates may be sent to the user terminal in 
3 ' N m • N a P • N sb bytes, where "3" is for the 24 total bits used to represent the I and Q 
components and N sb is the number of designated subbands. 

[0040] The user terminal also receives a downlink MIMO pilot transmitted by the 

access point (step 316) and derives an estimate of the downlink channel response, 
H dn (k) , for each of the designated subbands based on the received pilot (step 318). The 
user terminal then determines correction factors, K ap (£) and K m (k) , for each of the 
designated subbands based on the uplink and downlink channel response estimates, 
H up (£) and H dn (fc) (step 320). 

[0041] For the derivation of the correction factors, the downlink and uplink channel 

responses for each subband are assumed to be reciprocal, with gain/phase corrections to 
account for the differences in the transmit/receive chains at the access point and user 
terminal, as follows: 
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H up (*)K ut (&) = (H^K^C*)) 7 " , for k £ K , Eq (8) 

where K represents a set with all data subbands. Since only estimates of the effective 
downlink and uplink channel responses are available for the designated subbands during 
calibration, equation (8) may be rewritten as: 

Hup(*:)K ut (*) = (H dn (fc)K ap (£)) r , for k £ K' , Eq (9) 

where K' represents a set with all designated subbands. A correction vector k ut (£) 
may be defined to include only the N M diagonal elements of K ut (k) . Similarly, a 
correction vector k ap (£) may be defined to include only the N ap diagonal elements of 
K ap (fc). 

[0042] The correction factors K ap (£) and K ut (&) may be derived from the channel 

estimates H dn (£) and H up (fc) in various manners, including by a matrix-ratio 
computation and a minimum mean square error (MMSE) computation. Both of these 
computation methods are described in further detail below. Other computation methods 
may also be used, and this is within the scope of the invention. 

A. Matrix-Ratio Computation 
[0043] FIG. 4 is a flow diagram of an embodiment of a process 320a for deriving the 

correction vectors k ut (&) and k ap (&) from the downlink and uplink channel response 

estimates H up (&) and H^Cfc) using matrix-ratio computation. Process 320a may be 

used for step 320 in FIG. 3. 
[0044] Initially, an (N ut x N ap ) matrix C(£) is computed for each designated subband 

(step 412), as follows: 

~ T 

- {k) = W^ , forte r, Eq(10) 

where the ratio is taken element-by-element. Each element of C(&) may thus be 
computed as: 

c t j (k) = V jW , for i = {1 ... N ut } and j = {1 ... N ao } , Eq (1 1) 
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where h apiJ (k) and h ini J (k) are the >th (row, column) element of H„ p (£) and 
HanW , respectively, and c,^*) is the (i, j )-th element of C(£) . 
[0045] In an embodiment, the correction vector for the access point, k ap (k) , is defined 

to be equal to the mean of the normalized rows of C(k) and is derived by the steps in 
block 420. Each row of C(k) is first normalized by scaling each of the N ap elements in 
the row with the first element in the row (step 422). Thus, if 
c ; (fc) = [c t l (k) ... c liN (k)] is the z-th row of C(k) , then the normalized row £,(£) may 
be expressed as: 

c i (k) = [c il (k)/c iA (k) ... c fJ (*)/c u (*) ... c tiN Jk)/c ul (k)] . Eq(12) 
The mean of the normalized rows is then determined as the sum of the N ul normalized 
rows divided by N ut (step 424). The correction vector k ap (k) is set equal to this mean 
(step 426), which may be expressed as: 

L P (k) = ^-J^c i (k) , for keK'. Eq(13) 
N ut i=i 

Because of the normalization, the first element of k ap (k) is unity. 
[0046] In an embodiment, the correction vector for the user terminal, k ut (k) , is defined 

to be equal to the mean of the inverses of the normalized columns of C(k) and is 
derived by the steps in block 430. The j-th column of C(Jfc) is first normalized by 
scaling each element in the column with the j'-th element of the vector k ap (k) , which is 
denoted as K apJJ (k) (step 432). Thus, if c } (k) = [c hj (k) ... c„ mJ (k)f is the ;-th 
column of C(&) , then the normalized column c.j(k) may be expressed as: 

CjiV^jW/K^k) ... c Uj (k)/K a?JJ (k) ... c Nunj {k)IK &vJJ {k)f . Eq(14) 
The mean of the inverses of the normalized columns is then determined as the sum of 
the inverses of the N ap normalized columns divided by N ap (step 434). The correction 

vector k ut (k) is set equal to this mean (step 436), which may be expressed as: 
1 ^ 1 

kut(*) = — Z^T 77 ^orkeK, Eq(15) 
N ap Cj(k) 



WO 2004/039022 



PCT7US2003/034515 



13 

where the inversion of the normalized columns, Cj(k) , is performed element-wise. 

B. MMSE Computation 
[0047] For the MMSE computation, the correction factors K ap (k) and K ut (fc) are 

derived from the downlink and uplink channel response estimates and H up (&) 

such that the mean square error (MSE) between the calibrated downlink channel 
response and the calibrated uplink channel response is minimized. This condition may 
be expressed as: 

min |(H dn (£)K ap (*)) r -U up (k)K nt (k)\ 2 , for ke K, Eq (16) 

which may also be written as: 

min |K ap (k)El (k) - H up (*)K ut (*)| * , for k e K , 

* T 

where K ap (jk) = K ap (k) since K ap (£) is a diagonal matrix. 
[0048] Equation (16) is subject to the constraint that the lead element of K ap (£) is set 

equal to unity (i.e., K ap00 (fc) = 1). Without this constraint, the trivial solution would be 
obtained with all elements of the matrices K ap (k) and K ut (k) set equal to zero. In 

equation (16), a matrix Y(k) is first obtained as Y(A:) = K ap (yt)Hl(A:)-H up (A:)K ut (jt). 
The square of the absolute value is next obtained for each of the N ap • N ut entries of the 
matrix Y(k) . The mean square error (or the square error, since a divide by N ap • N ut is 
omitted) is then equal to the sum of all N ap • N ul squared values. 
[0049] The MMSE computation is performed for each designated subband to obtain the 

correction factors K ap (&) and K ut (&) for that subband. The MMSE computation for 
one subband is described below. For simplicity, the subband index, k, is omitted in the 
following description. Also for simplicity, the elements of the downlink channel 
response estimate H dn are denoted as {a tj } , the elements of the uplink channel response 

estimate H up are denoted as {b a } , the diagonal elements of the matrix K ap are denoted 
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as {m,.}, and the diagonal elements of the matrix K ut are denoted as {v,}, where 
/ = {1 ... iV flp }and; = {l ... NJ. 
[0050] The mean square error may be rewritten from equation (16), as follows: 

MSE = 2ZK".--Vvl 2 . Eq(17) 

again subject to the constraint «, = 1 . The minimum mean square error may be obtained 
by taking the partial derivatives of equation (17) with respect to u and v and setting the 
partial derivatives to zero. The results of these operations are the following equation 
sets: 

~ VP ' a *J = 0 ' for iG * 2 - N a P } . and Eq (18a) 

7=1 

-b ijVj )-b: =0 , for je {1 ... N ut } . Eq (18b) 

In equation (18a), u x =1 so there is no partial derivative for this case, and the index i 
runs from 2 through N ap . 
[0051] The set of (N ap +N ul -1) equations in equation sets (18a) and (18b) may be 

more conveniently expressed in matrix form, as follows: 

Ay = z , Eq(19) 

where 
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phjf o ... o 

° tW o ... 

0 ... 0 

- ° iK 



-b 0A cu, 



0 



- a 2N m b l Nm 



0 ... 0 



" U 2 " 




0 






0 


V, 


and z = 


0 


V 2 













[0052] The matrix A includes (Af flp + iV„, - 1) rows, with the first AT ap -l rows 

corresponding to the N ap -1 equations from equation set (18a) and the last N m rows 
corresponding to the N ut equations from equation set (18b). In particular, the first row 
of the matrix A is generated from equation set (18a) with i = 2, the second row is 
generated with i = 3 , and so on. The iV ap -th row of the matrix A is generated from 
equation set (18b) with j = 1, and so on, and the last row is generated with j = N ut . As 
shown above, the entries of the matrix A and the entries of the vector z may be 
obtained based on the entries in the matrices and H up . 

[0053] The correction factors are included in the vector y , which may be obtained as: 

y = A -1 z . Eq(20) 



WO 2004/039022 



PCT7US2003/034515 



[0054] The results of the MMSE computation are correction matrices K ap and K ut that 

minimize the mean square error in the calibrated downlink and uplink channel 
responses, as shown in equation (16). Since the matrices K ap and K ut are obtained 
based on the downlink and uplink channel response estimates, H dn and H up , the quality 
of the correction matrices K ap and K ut are thus dependent on the quality of the channel 
estimates and H up . The MIMO pilot may be averaged at the receiver to obtain 
more accurate estimates for and H up . 

[0055] The correction matrices, K ap and K ut , obtained based on the MMSE 

computation are generally better than the correction matrices obtained based on the 
matrix-ratio computation, especially when some of the channel gains are small and 
measurement noise can greatly degrade the channel gains. 

C. Post Computation 

[0056] Regardless of the particular computation method selected for use, after 

completion of the computation of the correction matrices, the user terminal sends to the 
access point the correction vectors for the access point, k ap (£), for all designated 
subbands. If 12-bit complex values are used for each correction factor in k ap (£) , then 
the correction vectors k ap (k) for all designated subbands may be sent to the access 
point in 3 ■ (N ap - 1) • N sb bytes, where "3" is for the 24 total bits used for the I and Q 
components and (N ap -1) results from the first element in each vector k ap (jfc) being 
equal to unity and thus not needing to be sent. If the first element is set to 
2 9 -1 = +511 , then 12 dB of headroom is available (since the maximum positive 12-bit 
signed value is 2" -1 = +2047 ), which would then allow gain mismatch of up to 12 dB 
between the downlink and uplink to be accommodated by 12-bit values. If the 
downlink and uplink match to within 12 dB and the first element is normalized to a 
value of 511, then the other elements should be no greater than 511-4 = 2044 in 
absolute value and can be represented with 12 bits. 
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[0057] A pair of correction vectors k ap (k) and k ut (k) is obtained for each designated 

subband. If the calibration is performed for fewer than all of the data subbands, then the 
correction factors for the "uncalibrated" subbands may be obtained by interpolating the 
correction factors obtained for the designated subbands. The interpolation may be 
performed by the access point to obtain the correction vectors k ap (k) , for k&K . 
Similarly, the interpolation may be performed by the user terminal to obtain the 
correction vectors k m (k) , for ke K . 

[0058] The access point and user terminal thereafter use their respective correction 

vectors k ap (&) and k ut (fc), or the corresponding correction matrices K ap (fc) and 

K ut (k) , for k G K , to scale modulation symbols prior to transmission over the wireless 
channel, as described below. The effective downlink channel that the user terminal sees 
would then be H cdn (k) = H dn (fc)K ap (*) . 

[0059] The calibration scheme described above, whereby a vector of correction factors 

is obtained for each of the access point and user terminal, allows "compatible" 
correction vectors to be derived for the access point when the calibration is performed 
by different user terminals. If the access point has already been calibrated (e.g., by one 
or more other user terminals), then the current correction vectors may be updated with 
the newly derived correction vectors. 

[0060] For example, if two user terminals simultaneously perform the calibration 

procedure, then the calibration results from these user terminals may be averaged to 
improve performance. However, calibration is typically performed for one user 
terminal at a time. So the second user terminal observes the downlink with the 
correction vector for the first user terminal already applied. In this case, the product of 
the second correction vector with the old correction vector may be used as the new 
correction vector, or a "weighted averaging" (described below) may also be used. The 
access point typically uses a single correction vector for all user terminals, and not 
different ones for different user terminals (although this may also be implemented). 
Updates from multiple user terminals or sequential updates from one user terminal may 
be treated in the same manner. The updated vectors may be directly applied (by a 
product operation). Alternatively, if some averaging is desired to reduce measurement 
noise, then weighted averaging may be used as described below. 
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[0061] Thus, if the access point uses correction vectors k apl (k) to transmit the MTMO 

pilot from which the user terminal determines new correction vectors k ap2 (k) , then the 
updated correction vectors k ap3 (fc) are the product of the current and new correction 
vectors. The correction vectors k apl (fe) and k ap2 (£) may be derived by the same or 
different user terminals. 

[0062] In one embodiment, the updated correction vectors are defined as 

k ap3 (&) = k apl (fc) -k ap2 (A:), where the multiplication is element-by-element. In another 
embodiment, the updated correction vectors may be redefined as 
ka P 3(£) = £. apl (^) -k ap2 (A:), where a is a factor used to provide weighted averaging (e.g., 
0 < a < 1 ). If the calibration updates are infrequent, then □ close to one might perform 
best. If the calibration updates are frequent but noisy, then a smaller value for □ is 
better. The updated correction vectors k ap3 (&) may then be used by the access point 
until they are updated again. 

[0063] As noted above, the calibration may be performed for fewer than all data 

subbands. For example, the calibration may be performed for every /i-th subband, 
where n may be determined by the expected response of the transmit/receive chains 
(e.g., n may be 2, 4, 8, 16, and so on). The calibration may also be performed for non- 
uniformly distributed subbands. For example, since there may be more filter roll-off at 
the edges of the passband, which may create more mismatch in the transmit/receive 
chains, more subbands near the band edges may be calibrated. In general, any number 
and any distribution of subbands may be calibrated, and this is within the scope of the 
invention. 

[0064] In the above description, the correction vectors k ap (£) and k ut (jfe) , for ke K', 

are derived by the user terminal, and the vectors k ap (k) are sent back to the access 
point. This scheme advantageously distributes the calibration processing among the 
user terminals for a multiple-access system. However, the correction vectors k ap (k) 
and k ut (k) may also be derived by the access point, which would then send the vectors 
k ut (k) back to the user terminal, and this is within the scope of the invention. 
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[0065] The calibration scheme described above allows each user terminal to calibrate 

its transmit/receive chains in real-time via over-the-air transmission. This allows user 
terminals with different frequency responses to achieve high performance without the 
need for tight frequency response specifications or to perform calibration at the factory. 
The access point may be calibrated by multiple user terminals to provide improved 
accuracy. 

D. Gain Considerations 

[0066] The calibration may be performed based on "normalized" gains for the downlink 

and uplink channels, which are gains given relative to the noise floor at the receiver. 
The use of the normalized gains allows the characteristics of one link (including the 
channel gains and SNR per eigenmode) to be obtained based on gain measurements for 
the other link, after the downlink and uplink have been calibrated. 

[0067] The access point and user terminal may initially balance their receiver input 

levels such that the noise levels on the receive paths for the access point and user 
terminal are approximately the same. The balancing may be done by estimating the 
noise floor, that is, finding a section of a received TDD frame (i.e., a unit of 
downlink/uplink transmission) that has a minimum average power over a particular time 
duration (e.g., one or two symbol periods). Generally, the time just before the start of 
each TDD frame is clear of transmissions, since any uplink data needs to be received by 
the access point and then a receive/transmit turnaround time is necessary before the 
access point transmits on the downlink. Depending on the interference environment, the 
noise floor may be determined based on a number of TDD frames. The downlink and 
uplink channel responses are then measured relative to this noise floor. More 
specifically, the channel gain for a given subband of a given transmit/receive antenna 
pair may first be obtained, for example, as the ratio of the received pilot symbol over the 
transmitted pilot symbol for that subband of that transmit/receive antenna pair. The 
normalized gain is then the measured gain divided by the noise floor. 

[0068] A large difference in the normalized gains for the access point and the 

normalized gains for the user terminal can result in the correction factors for the user 
terminal being greatly different from unity. The correction factors for the access point 
are close to unity because the first element of the matrix K ap is set to 1. 
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[0069] If the correction factors for the user terminal differ greatly from unity, then the 

user terminal may not be able to apply the computed correction factors. This is because 
the user terminal has a constraint on its maximum transmit power and may not be 
capable of increasing its transmit power for large correction factors. Moreover, a 
reduction in transmit power for small correction factors is generally not desirable, since 
this may reduce the achievable data rate. 

[0070] Thus, the user terminal can transmit using a scaled version of the computed 

correction factors. The scaled calibration factors may be obtained by scaling the 
computed correction factors by a particular scaling value, which may be set equal to a 
gain delta (difference or ratio) between the downlink and uplink channel responses. 
This gain delta can be computed as an average of the differences (or deltas) between the 
normalized gains for the downlink and uplink. The scaling value (or gain delta) used 
for the correction factors for the user terminal can be sent to the access point along with 
the computed correction factors for the access point. 

[0071] With the correction factors and the scaling value or gain delta, the downlink 

channel characteristics may be determined from the measured uplink channel response, 
and vice versa. If the noise floor at either the access point or the user terminal changes, 
then the gain delta can be updated, and the updated gain delta may be sent in a message 
to the other entity. 

[0072] In the above description, the calibration results in two sets (or vectors or 

matrices) of correction factors for each subband, with one set being used by the access 
point for downlink data transmission and the other set being used by the user terminal 
for uplink data transmission. The calibration may also be performed such that two sets 
of correction factors are provided for each subband, with one set being used by the 
access point for uplink data reception and the other set being used by the user terminal 
for downlink data reception. The calibration may also be performed such that one set 
of correction factors is obtained for each subband, and this set may be used at either the 
access point or the user terminal. In general, the calibration is performed such that the 
calibrated downlink and uplink channel responses are reciprocal, regardless of where 
correction factors are applied. 
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2. MIMO Pilot 

[0073] For the calibration, a MIMO pilot is transmitted on the uplink by the user 

terminal to allow the access point to estimate the uplink channel response, and a MIMO 
pilot is transmitted on the downlink by the access point to allow the user terminal to 
estimate the downlink channel response. The same or different MIMO pilots may be 
used for the downlink and uplink, and the MIMO pilots used are known at both the 
access point and user terminal. 

[0074] In an embodiment, the MIMO pilot comprises a specific OFDM symbol 

(denoted as "P") that is transmitted from each of the N T transmit antennas, where 
N T = N ap for the downlink and N T = N ul for the uplink. For each transmit antenna, the 
same P OFDM symbol is transmitted in each symbol period designated for MIMO pilot 
transmission. However, the P OFDM symbols for each antenna are covered with a 
different iV-chip Walsh sequence assigned to that antenna, where N>N ap for the 
downlink and N>N M for the uplink. The Walsh covering maintains orthogonality 
between the N T transmit antennas and allows the receiver to distinguish the individual 
transmit antennas. 

[0075] The P OFDM symbol includes one modulation symbol for each of the N sb 

designated subbands. The P OFDM symbol thus comprises a specific "word" of N sb 
modulation symbols that may be selected to facilitate channel estimation by the 
receiver. This word may also be defined to minimize the peak-to-average variation in 
the transmitted MIMO pilot. This may then reduce the amount of distortion and non- 
linearity generated by the transmit/receive chains, which may then result in improved 
accuracy for the channel estimation. 

[0076] For clarity, a specific MIMO pilot is described below for a specific MDVIO- 

OFDM system. For this system, the access point and user terminal each has four 
transmit/receive antennas. The system bandwidth is partitioned into 64 orthogonal 
subbands (i.e., N F =64), which are assigned indices of +31 to -32. Of these 64 
subbands, 48 subbands (e.g., with indices of ±{1, 6, 8, 20, 22, ... , 26}) are used 
for data, 4 subbands (e.g., with indices of +{7, 21}) are used for pilot and possibly 
signaling, the DC subband (with index of 0) is not used, and the remaining subbands are 
also not used and serve as guard subbands. This OFDM subband structure is described 
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in further detail in a document for IEEE Standard 802.11a and entitled "Part 11: 
Wireless LAN Medium Access Control (MAC) and Physical Layer (PHY) 
specifications: High-speed Physical Layer in the 5 GHz Band," September 1999, which 
is publicly available and incorporated herein by reference. 

[0077] The P OFDM symbol includes a set of 52 QPSK modulation symbols for the 48 

data subbands and 4 pilot subbands. This P OFDM symbol may be given as follows: 
P(real) = g • {0,0,0,0,0,0,-1,-1,-1,-1,1,1,1,-1,-1,1,-1,1,1,1,1,-1,-1,1,-1,1,-1,-1,-1,-1,1,-1, 

0,1,-1,-1,-1,-1,1,-1,-1,-1,-1,1,1,-1,-1,1,-1,-1,1,1,-1, l,-l,l,-l,l,-,0,0,0,0,0}, 
P(imag) = g ■ {0,0,0,0,0,0,-1,1,1,1,-1,-1,1,-1,1,1,1,-1,1,-1,-1,-1,-1,-1,-1,1,1,-1,1,1,-1,1, 

0,-1,-1,-1,-1,1,1,-1,1,-1,-1,1,-1,1,-1,1,1,1,-1,1,1,1,1,1,1,-1,-1,0,0,0,0,0}, 
where g is a gain for the pilot. The values within the { } bracket are given for subband 
indices -32 through -1 (for the first line) and 0 through +31 (for the second line). Thus, 
the first line for P(real) and P(imag) indicates that symbol (-1 - ;') is transmitted in 
subband -26, symbol (-1 + j) is transmitted in subband -25, and so on. The second line 
for P(real) and P(imag) indicates that symbol (1- ;') is transmitted in subband 1, 
symbol (-l-y) is transmitted in subband 2, and so on. Other OFDM symbols may 
also be used for the MTMO pilot. 

[0078] In an embodiment, the four transmit antennas are assigned Walsh sequences of 

W, = 1111 , W 2 = 1010, W 3 = 1100 , and W 4 = 1001 for the MEVIO pilot. For a given 
Walsh sequence, a value of "1" indicates that a P OFDM symbol is transmitted and a 
value of "0" indicates that a -P OFDM symbol is transmitted (i.e., each of the 52 
modulation symbols in P is inverted). 

[0079] Table 1 lists the OFDM symbols transmitted from each of the four transmit 

antennas for a MEVIO pilot transmission that spans four symbol periods. 



Table 1 



OFDM 
symbol 


Antenna 1 


Antenna 2 


Antenna 3 


Antenna 4 


1 


+P 


+P 


+P 


+P 


2 


+P 


-P 


+P 


-P 


3 


+P 


+P 


-P 


-P 


4 


+P 


-P 


-P 


+P 
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For longer MIMO pilot transmission, the Walsh sequence for each transmit antenna is 
simply repeated. For this set of Walsh sequences, the MIMO pilot transmission occurs 
in integer multiples of four symbol periods to ensure orthogonality among the four 
transmit antennas. 

[0080] The receiver may derive an estimate of the channel response based on the 

received MIMO pilot by performing the complementary processing. In particular, to 
recover the pilot sent from transmit antenna i and received by receive antenna j, the pilot 
received by receive antenna j is first processed with the Walsh sequence assigned to 
transmit antenna i in a complementary manner to the Walsh covering performed at the 
transmitter. The decovered OFDM symbols for all N ps symbol periods for the MDVIO 
pilot are then accumulated, where the accumulation is performed individually for each 
of the 52 subbands used to carry the MIMO pilot. The results of the accumulation is 
h u (k) , for k = ± {1, 26} , which is an estimate of the effective channel response from 
transmit antenna i to receive antenna j (i.e., including the responses for the 
transmit/receive chains) for the 52 data and pilot subbands. 

[0081] The same processing may be performed to recover the pilot from each transmit 

antenna at each receive antenna. The pilot processing provides N ap • N M values that are 

the elements of the effective channel response estimate, H up (A;) or H dn (£) , for each of 
the 52 subbands. 

[0082] The channel estimation described above may be performed by both the access 

point and the user terminal during calibration to obtain the effective uplink channel 
response estimate H up (&) and the effective downlink channel response estimate 

H dn (&) , respectively, which are then used to derive the correction factors as described 
above. 

3. Spatial Processing 

[0083] The correlation between the downlink and uplink channel responses may be 

exploited to simplify the channel estimation and spatial processing at the access point 
and user terminal for TDD MIMO and MIMO-OFDM systems. This simplification is 
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possible after calibration has been performed to account for differences in the 
transmit/receive chains. As noted above, the calibrated channel responses are: 

H cdn (k) = H (In (&)K ap (k) , for the downlink, and Eq (21a) 

S cup (*) = (H^ (fc)K ap (k)) T = H up (fc)K ut (k) , for the uplink. Eq (21b) 

The last equality in equation (21b) comes from using the relationship between the 
effective downlink and uplink channel responses, H up (k) = (K^ (£)H dn (£)K ap (Jc)) T . 

[0084] The channel response matrix H(£) for each subband may be "diagonalized" to 

obtain the Ns eigenmodes for that subband. This may be achieved by performing either 
singular value decomposition on the channel response matrix H(£) or eigenvalue 
decomposition on the correlation matrix of H(fc) , which is R(k) = H H (fc)H(Jfc) . 

[0085] The singular value decomposition of the calibrated uplink channel response 

matrix, H cup (£) , may be expressed as: 

= U ap (*)2(fc)V u *; (*) , for ke K , Eq (22) 

where U ap (fc) is an (N ut x N ut ) unitary matrix of left eigenvectors of H cup (k) ; 
S(£) is an (N ut xN ap ) diagonal matrix of singular values of H cup (&) ; and 
V ut (fc) is an (N ap xN ap ) unitary matrix of right eigenvectors of H cup (£) . 
[0086] A unitary matrix M is characterized by the property M H M = I , where I is the 

identity matrix. Correspondingly, the singular value decomposition of the calibrated 
downlink channel response matrix, H cdn (k) , may be expressed as: 

H cdn W = V; t WZWUl p W , for ke K. Eq (23) 

The matrices V* (k) and lf ap (k) are thus also matrices of left and right eigenvectors, 
respectively, of H cdn (£) . The matrices V m (Jk), V* t (Jfc), Y T ul (k), and \"(k) are 
different forms of the matrix V ut (k) , and the matrices U ap (k) , U* p (k) , U ap (k) , and 
Uf p (*) are also different forms of the matrix U ap (£) . For simplicity, reference to the 
matrices U ap (&) and V ut (&) in the following description may also refer to their various 
other forms. The matrices U ap (A;) and V ut (A:) are used by the access point and user 
terminal, respectively, for spatial processing and are denoted as such by their subscripts. 
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[0087] The singular value decomposition is described in further detail by Gilbert Strang 

entitled "Linear Algebra and Its Applications," Second Edition, Academic Press, 1980. 

[0088] The user terminal can estimate the calibrated downlink channel response based 

on a MIMO pilot sent by the access point. The user terminal may then perform the 
singular value decomposition of the calibrated downlink channel response estimate 
H cdn (&), for ke K , to obtain the diagonal matrices t(k) and the matrices V* t (£) of 
left eigenvectors of H cdn (£). This singular value decomposition may be given as 

HcdnW = Y ut (fc)2(fc)ulp(&) , where the hat (" A ") above each matrix indicates that it is 
an estimate of the actual matrix. 
[0089] Similarly, the access point can estimate the calibrated uplink channel response 

based on a MIMO pilot sent by the user terminal. The access point may then perform 
the singular value decomposition of the calibrated uplink channel response estimate 

H cup (A:) , for k E K , to obtain the diagonal matrices l(k) and the matrices U ap (£) of 
left eigenvectors of H cup (£), for ke K . This singular value decomposition may be 
given as H cup (£) = U ap (fc)i(*)V u "(fc). 
[0090] Because of the reciprocal channel and the calibration, the singular value 

decomposition only need to be performed by either the user terminal or the access point 
to obtain both matrices V ut (fc) and U ap (fc) . If performed by the user terminal, then the 
matrices V ut (k) are used for spatial processing at the user terminal and the matrices 
U ap (&) may be sent back to the access point. 

[0091] The access point may also be able to obtain the matrices U ap (&) and l(k) based 

on a steered reference sent by the user terminal. Similarly, the user terminal may also 
be able to obtain the matrices V m (£) and t(k) based on a steered reference sent by the 
access point. The steered reference is described in detail in the aforementioned 
provisional U.S. Patent Application Serial No. 60/421,309. 

[0092] The matrices U ap (£) and l(k) may be used to transmit independent data 

streams on the N s eigenmodes of the MIMO channel, where N s < min . The 
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spatial processing to transmit multiple data streams on the downlink and uplink is 
described below. 

A. Uplink Spatial Processing 

[0093] The spatial processing by the user terminal for an uplink transmission may be 

expressed as: 

x up (^) = K ut (fe)V ut (^)s up (fc) JorkeK, Eq (24) 

where x up (k) the transmit vector for the uplink for the k-th subband; and 

s up (£) is a "data" vector with up to N s non-zero entries for the modulation 
symbols to be transmitted on the N s eigenmodes of the k-th subband. 
[0094] Additional processing may also be performed on the modulation symbols prior 

to transmission. For example, channel inversion may be applied across the data 
subbands (e.g., for each eigenmode) such that the received SNR is approximately equal 
for all data subbands. The spatial processing may then be expressed as: 

x up (*) = K ut (k)V ut (kW ap (k)s» p (k) , for k e K , Eq (25) 

where W up (k) is a matrix with weights for the (optional) uplink channel inversion. 

[0095] The channel inversion may also be performed by assigning transmit power to 

each subband before the modulation takes place, in which case the vector s up (fc) 
includes the channel inversion coefficients and the matrix W up (k) can be omitted from 
equation (25). In the following description, the use of the matrix W up (£) in an 
equation indicates that the channel inversion coefficients are not incorporated into the 
vector s up (&) . The lack of the matrix W up (&) in an equation can indicate either (1) 
channel inversion is not performed or (2) channel inversion is performed and 
incorporated into the vector s up (£) . 

[0096] Channel inversion may be performed as described in the aforementioned 

provisional U.S. Patent Application Serial No. 60/421,309 and in U.S. Patent 
Application Serial No. 10/229,209, entitled "Coded MEMO Systems with Selective 
Channel Inversion Applied Per Eigenmode," filed August 27, 2002, assigned to the 
assignee of the present application and incorporated herein by reference. 
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[0097] The received uplink transmission at the access point may be expressed as: 

r up (&) =U np (k)x up (k) + n(k) JovkeK, Eq(26) 

where r up (k) is the received vector for the uplink for the k-th subband; 

n(k) is additive white Gaussian noise (AWGN) for the k-th subband; and 
x up (&) is as shown in equation (24). 
[0098] The spatial processing (or matched filtering) at the access point for the received 

uplink transmission may be expressed as: 

M*) =t~\k)tll(k)r up (k) 

= t~\k)t^ ? (k)^(k)±(k)s^(k) + n(k)) JovkeK, Eq (27) 
= s up (fc) + n(fc) 

where s up (fc) is an estimate of the vector s up (£) transmitted by the user terminal on the 
uplink, and n(k) is the post-processed noise. Equation (27) assumes that channel 
inversion was not performed at the transmitter and that the received vector r up (k) is as 
shown in equation (26). 

B. Downlink Spatial Processing 
[0099] The spatial processing by the access point for a downlink transmission may be 

expressed as: 

XdnW=Ka P Wul p (fe)s dn (A:) JovkeK, Eq (28) 

where x dn (^) is the transmit vector and s in (k) is the data vector for the downlink. 
[00100] Again, additional processing (e.g., channel inversion) may also be performed on 

the modulation symbols prior to transmission. The spatial processing may then be 
expressed as: 

x d „ (*) = K ap (k)tf ap (£) W dn (*)s dn (k) JovkeK, Eq (29) 

where W dn (k) is a matrix with weights for the (optional) downlink channel inversion. 
[00101] The received downlink transmission at the user terminal may be expressed as: 
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r dn (k) = H, n (k)x ta (k) + n(k) , for k £ K , Eq (30) 

where (A:) is the transmit vector as shown in equation (28). 
[00102] The spatial processing (or matched filtering) at the user terminal for the received 

downlink transmission may be expressed as: 

§„„(*) =l" 1 w£(^)r dn (fc) 

= t~\k)VM(fAk)tm aa (k)+n(k)) JorkeK. Eq(31) 
= s dn (fe) + n(fc) 

Equation (31) assumes that channel inversion was not performed at the transmitter and 
that the received vector (k) is as shown in equation (30). 
[00103] Table 2 summarizes the spatial processing at the access point and user terminal 

for data transmission and reception. Table 2 assumes that the additional processing by 
W(£) is performed at the transmitter. However, if this additional processing is not 
performed, then W(k) can be regarded as the identify matrix. 



Table 2 





Uplink 


Downlink 


User 
Terminal 


Transmit : 
x up (*) = K ut (fc)V ut (^)W up (A:)s up (^) 


Receive : 
s ia (k) = t~\k)V ut (k)r tiD (k) j 


Access 
Point 


Receive : 
| up (fc) = r 1 (fc)Ur p (^)r up W 


Transmit : 
XdnW = K ap (A:)U ap (^)W dn (A:)s dn (fc) 



[00104] In the above description and as shown in Table 2, the correction matrices 

K ap (k) and K ut (k) are used for the transmit spatial processing at the access point and 
user terminal, respectively. This can simplify the overall spatial processing since the 
modulation symbols may need to be scaled anyway (e.g., for channel inversion) and the 
correction matrices K ap (&) and K nt (fc) may be combined with the weight matrices 
W dn (£) and W up (fc) to obtain gain matrices G dn (&) and G ap (k) , where 
GdnW = W dn (£)K ap (£) and G up (*) = W up (Jfc)K tt (k) . The processing may also be 



WO 2004/039022 



PCT7US2003/034515 



29 

performed such that the correction matrices are used for the receive spatial processing 
(instead of the transmit spatial processing). 

4. MIMO-OFDM System 

[00105] FIG. 5 is a block diagram of an embodiment of an access point 502 and a user 

terminal 504 within a TDD MIMO-OFDM system. For simplicity, the following 
description assumes that the access point and user terminal are each equipped with four 
transmit/receive antennas. 

[00106] On the downlink, at access point 502, a transmit (TX) data processor 510 

receives traffic data (i.e., information bits) from a data source 508 and signaling and 
other information from a controller 530. TX data processor 510 formats, codes, 
interleaves, and modulates (i.e., symbol maps) the data to provide a stream of 
modulation symbols for each eigenmode used for data transmission. A TX spatial 
processor 520 receives the modulation symbol streams from TX data processor 510 and 
performs spatial processing to provide four streams of transmit symbols, one stream for 
each antenna. TX spatial processor 520 also multiplexes in pilot symbols as appropriate 
(e.g., for calibration). 

[00107] Each modulator (MOD) 522 receives and processes a respective transmit symbol 

stream to provide a corresponding stream of OFDM symbols. Each OFDM symbol 
stream is further processed by a transmit chain within modulator 522 to provide a 
corresponding downlink modulated signal. The four downlink modulated signals from 
modulator 522a through 522d are then transmitted from four antennas 524a through 
524d, respectively. 

[00108] At user terminal 504, antennas 552 receive the transmitted downlink modulated 

signals, and each antenna provides a received signal to a respective demodulator 
(DEMOD) 554. Each demodulator 554 (which includes a receive chain) performs 
processing complementary to that performed at modulator 522 and provides received 
symbols. A receive (RX) spatial processor 560 then performs spatial processing on the 
received symbols from all demodulators 554 to provide recovered symbols, which are 
estimates of the modulation symbols sent by the access point. During calibration, RX 
spatial processor 560 provides a calibrated downlink channel estimate, H cdn (£), based 
on the MFMO pilot transmitted by the access point. 
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[00109] An RX data processor 570 processes (e.g., symbol demaps, deinterleaves, and 

decodes) the recovered symbols to provide decoded data. The decoded data may 
include recovered traffic data, signaling, and so on, which are provided to a data sink 
572 for storage and/or a controller 580 for further processing. During calibration, RX 
data processor 570 provides the calibrated uplink channel estimate, H cup (&), which is 
derived by the access point and sent on the downlink. 

[00110] Controllers 530 and 580 control the operation of various processing units at the 

access point and user terminal, respectively. During calibration, controller 580 may 
receive the channel response estimates H cdn (£) and H cup (&) , derive the correction 
matrices K ap (£) and K ut (fc), provide the matrices K ut (&) to a TX spatial processor 
592 for uplink transmission, and provide the matrices K ap (k) to a TX data processor 
590 for transmission back to the access point. Memory units 532 and 582 store data and 
program codes used by controllers 530 and 580, respectively. 

[00111] The processing for the uplink may be the same or different from the processing 
for the downlink. Data and signaling are processed (e.g., coded, interleaved, and 
modulated) by a TX data processor 590 and further spatially processed by TX spatial 
processor 592, which multiplexes in pilot symbols during calibration. The pilot and 
modulation symbols are further processed by modulators 554 to generate uplink 
modulated signals, which are then transmitted via antennas 552 to the access point. 

[00112] At access point 1 10, the uplink modulated signals are received by antennas 524, 

demodulated by demodulators 522, and processed by an RX spatial processor 540 and 
an RX data processor 542 in a complementary to that performed at the user terminal. 
During calibration, RX spatial processor 560 also provides a calibrated uplink channel 
estimate, H cup (k) , based on the MIMO pilot transmitted by the user terminal. The 
matrices H cup (k) are received by controller 530 and then provided to TX data processor 
510 for transmission back to the user terminal. 

[00113] FIG. 6 is a block diagram of a TX spatial processor 520a, which may be used 

for TX spatial processors 520 and 592 in FIG. 5. For simplicity, the following 
description assumes that all four eigenmodes are selected for use. 
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[00114] Within processor 520a, a demultiplexer 632 receives four modulation symbol 
steams (denoted as through s 4 (n)) to be transmitted on four eigenmodes, 

demultiplexes each stream into N D substreams for the N D data subbands, and provides 
four modulation symbol substreams for each data subband to a respective TX subband 
spatial processor 640. Each processor 640 performs the processing shown in equation 
(24), (25), (28), or (29) for one subband. 

[00115] Within each TX subband spatial processor 640, the four modulation symbol 

substreams (denoted as s^k) through s 4 (k)) are provided to four multipliers 642a 
through 642d, which also receive the gains g,(fc), g 2 (k), g 3 (k), and g 4 (&) for the four 
eigenmodes of the associated subband. For the downlink, the four gains for each data 
subband are the diagonal elements of the corresponding matrix G dn (£), where 

G dn (fc) = K ap (&) or G dn (£) = W dn (fc)K ap (£) . For the uplink, the gains are the diagonal 
elements of the matrix G up (&), where G up (k) = K ut (k) or G np (k) = W up (£)K ut (fc) . 
Each multiplier 642 scales its modulation symbols with its gain g m (k) to provide scaled 
modulation symbols. Multipliers 642a through 642d provides four scaled modulation 
symbol substreams to four beam-formers 650a through 650d, respectively. 
[00116] Each beam-former 650 performs beam-forming to transmit one symbol 

substream on one eigenmode of one subband. Each beam-former 650 receives one 
scaled symbol substream s m (k) and performs beam-forming using the eigenvector 
v m (k) for the associated eigenmode. Within each beam-former 650, the scaled 
modulation symbols are provided to four multipliers 652a through 652d, which also 
receive four elements, v ffll (A:), v m2 (k), v m3 (k), and v mi (k), of the eigenvector \ m {k) 
for the associated eigenmode. The eigenvector v m (k) is the m-th column of the matrix 
U ap (£) for the downlink and is the m-th column of the matrix V ut (fe) for the uplink. 
Each multiplier 652 then multiplies the scaled modulation symbols with its eigenvector 
value v m j (k) to provide "beam-formed" symbols. Multipliers 652a through 652d 
provide four beam-formed symbol substreams (which are to be transmitted from four 
antennas) to summers 660a through 660d, respectively. 
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[00117] Each summer 660 receives and sums four beam-formed symbols for the four 

eigenmodes for each symbol period to provide a preconditioned symbol for an 
associated transmit antenna. Summers 660a through 660d provides four substreams of 
preconditioned symbols for four transmit antennas to buffers/multiplexers 670a through 
670d, respectively. 

[00118] Each buffer/multiplexer 670 receives pilot symbols and the preconditioned 

symbols from TX subband spatial processors 640 for the N D data subbands. Each 
buffer/multiplexer 670 then multiplexes pilot symbols, preconditioned symbols, and 
zeros for the pilot subbands, data subbands, and unused subbands, respectively, to form 
a sequence of N F transmit symbols for that symbol period. During calibration, pilot 
symbols are transmitted on the designated subbands. Multipliers 668a through 668d 
cover the pilot symbols for the four antennas with Walsh sequences W, through W 4 , 
respectively, assigned to the four antennas, as described above and shown in Table 1. 
Each buffer/multiplexer 670 provides a stream of transmit symbols x t (ri) for one 
transmit antenna, where the transmit symbol stream comprises concatenated sequences 
of N F transmit symbols. 

[00119] The spatial processing and OFDM modulation is described in further detail in 

the aforementioned provisional U.S. Patent Application Serial No. 60/421,309. 

[00120] In various embodiments of the invention as described herein, peer-peer 

communication between the various user terminals (UTs or STAs) in the same basic 
service set (BSS) or different BSSs can be implemented as described below. The UTs 
or STAs that calibrate with a single access point (AP) are members of a basic service set 
(BSS). The single access point is a common node to all UTs in the BSS. The calibration 
methods as described above facilitate the following types of communication: 

(i) A UT in the BSS can use TX steering to communicate directly with the AP 
on the uplink (UL) and the AP can use TX steering to communicate with the UTs on the 
downlink (DL). 

(ii) A UT in the BSS can communicate directly with another UT in the same 
BSS using steering. In this case, this peer-peer communication has to be bootstrapped 
because neither UT knows the channel between them. In various embodiments, the 
bootstrap procedure works as follows: 
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- The initiator of the peer-peer link is the designate AP (DAP), and the other UT 
is the designated UT (DUT). 

- The DAP sends MIMO pilot to the DUT along with a request to establish link, 
which contains the BSS ID plus the DAP ID. The request needs to be sent in a 
common mode (i.e. Tx diversity). 

- The DUT responds by sending back steered MIMO pilot plus an 
acknowledgement which contains the DUT ID, its BSS ID, and some rate 
indicator for the DAP to use. 

- The DAP can then use steering on the DL and the DUT can use steering on the 
UL. Rate control and tracking can be accommodated by breaking the 
transmissions into DL and UL segments with sufficient time between them to 
allow for processing. 

(iii) UTs that belong to one BSS (e.g., BSS1) can steer to UTs that belong to 
another BSS (e.g., BSS2), even though each has calibrated with a different AP. 
However, there will be a phase rotation ambiguity (per subband) in this case. This is 
because the calibration procedure as described above establishes a reference which is 
unique to the AP it has calibrated with. The reference is a complex constant, 

a(k,j)= 8AFnc( - 0) 
W(0) 

where k is the subband index and j is the AP index and 0 is the index of the reference 
antenna (e.g., antenna 0) used on the AP. In one embodiment, this constant is common 
to all UTs in a given BSS, but is independent for different BSSs. 

[00121] As a result, when a UT from BSS1 communicates with a UT in BSS2, steering 

without correction or compensation for this constant may result in a phase rotation and 
amplitude scaling of the entire eigensystem. The phase rotation can be determined 
through the use of pilot (steered and unsteered) and removed in the receivers of each 
respective UT. In one embodiment, the amplitude correction or compensation can 
simply be an SNR scaling and can be removed by estimation of the noise floor at each 
receiver, which may impact rate selection. 

[00122] In various embodiments, the peer-peer exchange between UTs that belong to 

different BSSs may work as follows: 

- The initiator of the peer-peer link (e.g., UT in BSS1) is the designate AP 
(DAP), and the other UT (e.g., UT in BSS2) is the designated UT (DUT). 
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- The DAP sends MIMO pilot to the DUT along with a request to establish link, 
which contains the respective BSS ID plus the DAP ID. The request needs to be sent in 
a common mode (i.e. Tx diversity). 

- The DUT responds by sending back steered MIMO pilot plus an 
acknowledgement which contains DUT ID, its BSS ID, and some rate indicator for the 
DAP to use. 

- The DAP receiver (Rx) can estimate the phase rotation on the uplink (UL) and 
apply the correction constant to each subband. The DAP can then use steering on the 
downlink (DL) but needs to include a preamble of steered reference on at least the first 
steered packet to allow the DUT receiver (Rx) to correct or compensate for the phase 
rotation on the DL for each subband. Subsequent DL transmissions may not require a 
steered reference preamble. Rate control and tracking can be accommodated by 
breaking the transmissions into DL and UL segments with sufficient time between them 
to allow for processing. 

[00123] The calibration techniques described herein may be implemented by various 
means. For example, these techniques may be implemented in hardware, software, or a 
combination thereof. For a hardware implementation, the techniques may be 
implemented at the access point and user terminal within one or more application 
specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal 
processing devices (DSPDs), programmable logic devices (PLDs), field programmable 
gate arrays (FPGAs), processors, controllers, micro-controllers, microprocessors, other 
electronic units designed to perform the functions described herein, or a combination 
thereof. 

[00124] For a software implementation, the calibration techniques may be implemented 

with modules (e.g., procedures, functions, and so on) that perform the functions 
described herein. The software codes may be stored in a memory unit (e.g., memory 
units 532 and 582 in FIG. 5) and executed by a processor (e.g., controllers 530 and 580, 
as appropriate). The memory unit may be implemented within the processor or external 
to the processor, in which case it can be communicatively coupled to the processor via 
various means as is known in the art. 

[00125] Headings are included herein for reference and to aid in locating certain 

sections. These headings are not intended to limit the scope of the concepts described 
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therein under, and these concepts may have applicability in other sections throughout 
the entire specification. 

[00126] The previous description of the disclosed embodiments is provided to enable any 

person skilled in the art to make or use the present invention. Various modifications to 
these embodiments will be readily apparent to those skilled in the art, and the generic 
principles defined herein may be applied to other embodiments without departing from 
the spirit or scope of the invention. Thus, the present invention is not intended to be 
limited to the embodiments shown herein but is to be accorded the widest scope 
consistent with the principles and novel features disclosed herein. 



WHAT IS CLAIMED IS: 
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CLAIMS 

1. A method for calibrating downlink and uplink channels in a wireless 
communication system, comprising: 

obtaining an estimate of a downlink channel response; 
obtaining an estimate of an uplink channel response; 

determining first and second sets of correction factors based on the estimates of 
the downlink and uplink channel responses; and 

calibrating the downlink channel and uplink channel based on the first and 
second sets of correction factors, respectively, to form a calibrated downlink channel 
and a calibrated uplink channel. 

2. The method of claim 1, wherein the first set of correction factors is used 
to scale symbols prior to transmission on the downlink channel and the second set of 
correction factors is used to scale symbols prior to transmission on the uplink channel. 

3. The method of claim 1, wherein the first set of correction factors is used 
to scale symbols received on the downlink channel and the second set of correction 
factors is used to scale symbols received on the uplink channel. 

4. The method of claim 1, wherein the first and second sets of correction 
factors are determined based on the following equation: 

H up K ut =(H dn K ap ) r , 
where is a matrix for the estimate of the downlink channel response, 
H up is a matrix for the estimate of the uplink channel response, 
K ap is a matrix for the first set of correction factors, 

K ut is a matrix for the second set of correction factors, and 
" T " denotes a transpose. 

5. The method of claim 4, wherein determining the first and second sets of 
correction factors includes: 
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computing a matrix C as an element-wise ratio of the matrix H up over the 
matrix , and 

deriving the matrices K ap and K ut based on the matrix C . 

6. The method of claim 5, wherein the deriving the matrix K ut includes 
normalizing each of a plurality of rows of the matrix C , and 

determining a mean of the plurality of normalized rows of the matrix C , and 
wherein the matrix K ut is formed based on the mean of the plurality of normalized 
rows. 

7. The method of claim 5, wherein the deriving the matrix K ap includes 
normalizing each of a plurality of columns of the matrix C , and 
determining a mean of inverses of the plurality of normalized columns of the 

matrix C , and wherein the matrix K ap is formed based on the mean of the inverses of 
the plurality of normalized columns. 

8. The method of claim 4, wherein the matrices K ut and K ap are derived 
based on a minimum mean square error (MMSE) computation. 

9. The method of claim 8, wherein the MMSE computation minimizes a 
mean square error (MSE) given as 

|H up K ut -(H dn K ap )T- 

10. The method of claim 1, further comprising: 

determining a scaling value indicative of an average difference between the 
estimate of the downlink channel response and the estimate of the uplink channel 
response. 
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11. The method of claim 1, wherein the estimates for the downlink and 
uplink channel responses are normalized to account for receiver noise floor. 

12. The method of claim 1, wherein the determining is performed at a user 
terminal. 



13. The method of claim 4, wherein a first set of matrices of correction 
factors for the downlink channel is determined for a first set of subbands, the method 
further comprising: 

interpolating the first set of matrices to obtain a second set of matrices of 
correction factors for the downlink channel for a second set of subbands. 

14. The method of claim 1, wherein the estimates of the downlink and uplink 
channel responses are each obtained based on a pilot transmitted from a plurality of 
antennas and orthogonalized with a plurality of orthogonal sequences. 

15. The method of claim 1 wherein the estimate of the uplink channel, 
response is obtained based on a pilot transmitted on the uplink channel and wherein the 
estimate of the downlink channel response is obtained based on a pilot transmitted on 
the downlink channel. 



16. The method of claim 1, wherein the TDD system is a multiple-input 
multiple-output (MIMO) system. 

17. The method of claim 1, wherein the TDD system utilizes orthogonal 
frequency division multiplexing (OFDM). 

18. A method for calibrating downlink and uplink channels in a wireless 
time division duplexed (TDD) multiple-input multiple-output (MIMO) communication 
system, comprising: 

transmitting a pilot on the uplink channel; 

obtaining an estimate of an uplink channel response derived based on the pilot 
transmitted on the uplink channel; 
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receiving a pilot on the downlink channel; 

obtaining an estimate of a downlink channel response derived based on the pilot 
received on the downlink channel; and 

determining first and second sets of correction factors based on the estimates of 
the downlink and uplink channel responses, wherein a calibrated downlink channel is 
formed by using the first set of correction factors for the downlink channel and a 
calibrated uplink channel is formed by using the second set of correction factors for the 
uplink channel. 

19. The method of claim 18, wherein the first and second sets of correction 
factors are determined based on a minimum mean square error (MMSE) computation. 

20. The method of claim 18, wherein the first and second sets of correction 
factors are determined based on a matrix-ratio computation. 

21. The method of claim 18, wherein the first set of correction factors is 
updated based on calibration with a plurality of user terminals. 

22. The method of claim 18, further comprising: 

scaling symbols with the first set of correction factors prior to transmission on 
the downlink. 

23. The method of claim 1 8, further comprising: 

scaling symbols with the second set of correction factors prior to transmission on 
the uplink channel. 

24. An apparatus in a wireless time division duplexed (TDD) multiple-input 
multiple-output (MIMO) communication system, comprising: 

means for obtaining an estimate of a response of a downlink channel; 
means for obtaining an estimate of a response of an uplink channel; and 
means for determining first and second sets of correction factors based on the 
estimates of the downlink and uplink channel responses, wherein a calibrated downlink 
channel is formed by using the first set of correction factors for the downlink channel 
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and a calibrated uplink channel is formed by using the second set of correction factors 
for the uplink channel. 

25. A user terminal in a wireless time division duplexed (TDD) 
communication system, comprising: 

an TX spatial processor operative to transmit a first pilot on an uplink channel; 

an RX spatial processor operative to receive a second pilot on a downlink 
channel and derive an estimate of a downlink channel response based on the received 
second pilot, and to receive an estimate of an uplink channel response derived based on 
the transmitted first pilot; and 

a controller operative to determine first and second sets of correction factors 
based on the estimates of the downlink and uplink channel responses, wherein a 
calibrated downlink channel is formed by using the first set of correction factors for the 
downlink channel and a calibrated uplink channel is formed by using the second set of 
correction factors for the uplink channel 

26. The user terminal of claim 25, wherein the controller is further operative 
to determine the first and second sets of correction factor based on a minimum mean 
square error (MMSE) computation. 

27. The user terminal of claim 25, wherein the controller is further operative 
to determine the first and second sets of correction factor based on a matrix-ratio 
computation. 

28. A method for communication in a wireless system, comprising: 
calibrating one or more communication links between a plurality of user stations 

and one or more access points, based on one or more sets of correction factors derived 
from estimates of channel responses associated with the one or more communication 
links, the plurality of user stations including a first user station and a second user 
station; and 

establishing communication between the first and second user stations using 
steering without performing calibration between the first and second user stations. 
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29. The method of claim 28 wherein establishing the communication 
between the first and second user stations comprises: 

sending, from the first user station, a pilot and a request to establish a 
communication link with the second user station; 

sending, from the second user station, a steered pilot and an acknowledgment in 
response to receiving the pilot and the request from first user station; 

transmitting information between the first and second user stations using 
steering based on the steered pilot. 

30. The method of claim 29 wherein the request to establish the 
communication comprises an identifier of a basic service set to which the first user 
station belongs and an identifier of the first user station. 

31. The method of claim 29 wherein the acknowledgment comprises an 
identifier of the second user station, an identifier of a basic service set to which the 
second user station belongs, and a data rate indicator. 

32. The method of claim 28 wherein the one or more access points includes a 
first access point associated with a first basic service set (BSS) and a second access 
point associated with a second BSS, wherein the first user station is calibrated with 
respect to the first access point and the second user station is calibrated with respect to 
the second access point, and wherein establishing the communication between the first 
and second user stations comprises: 

sending, from the first user station, a pilot and a request to establish a 
communication link with the second user station; 

sending, from the second user station, a steered pilot and an acknowledgment in 
response to receiving the pilot and the request from first user station; and 

transmitting information between the first and second user stations using 
steering that is adjusted to compensate for a phase rotation caused by calibration of the 
first and second user stations with respect to different access points. 

33. The method of claim 32 wherein the phase rotation is determined based 
on the steered pilot received from the second user station. 
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34. An apparatus for communication in a wireless system, comprising: 
means for alibrating one or more communication links between a plurality of 

user stations and one or more access points, based on one or more sets of correction 
factors derived from estimates of channel responses associated with the one or more 
communication links, the plurality of user stations including a first user station and a 
second user station; and 

means for establishing communication between the first and second user stations 
using steering without performing calibration between the first and second user stations. 

35. The apparatus of claim 34 wherein establishing the communication 
between the first and second user stations comprises: 

means for sending, from the first user station, a pilot and a request to establish a 
communication link with the second user station; 

means for sending, from the second user station, a steered pilot and an 
acknowledgment in response to receiving the pilot and the request from first user 
station; 

means for transmitting information between the first and second user stations 
using steering based on the steered pilot. 

36. The apparatus of claim 35 wherein the request to establish the 
communication comprises an identifier of a basic service set to which the first user 
station belongs and an identifier of the first user station. 

37. The apparatus of claim 35 wherein the acknowledgment comprises an 
identifier of the second user station, an identifier of a basic service set to which the 
second user station belongs, and a data rate indicator. 

38. The apparatus of claim 34 wherein the one or more access points 
includes a first access point associated with a first basic service set (BSS) and a second 
access point associated with a second BSS, wherein the first user station is calibrated 
with respect to the first access point and the second user station is calibrated with 
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respect to the second access point, and wherein establishing the communication between 
the first and second user stations comprises: 

sending, from the first user station, a pilot and a request to establish a 
communication link with the second user station; 

sending, from the second user station, a steered pilot and an acknowledgment in 
response to receiving the pilot and the request from first user station; and 

transmitting information between the first and second user stations using 
steering that is adjusted to compensate for a phase rotation caused by calibration of the 
first and second user stations with respect to different access points. 

39. The apparatus of claim 38 wherein the phase rotation is determined 
based on the steered pilot received from the second user station. 
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